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Abstract 

The European XFEL, currently under construction, will produce a coherent X-ray pulse every 222 ns in pulse 
trains of up to 2700 pulses. In conjunction with the fast 2D area detectors currently under development, it 
will be possible to perform X-ray Photon Correlation Spectroscopy (XPCS) experiments on sub-microsecond 
timescales with non-ergodic systems. 

A case study for the Adaptive Gain Integrating Pixel Detector (AGIPD) at the European XFEL employ- 
ing the intensity autocorrelation technique was performed using the detector simulation tool HORUS. As 
optimum results from XPCS experiments are obtained when the pixel size approximates the (small) speckle 
size, the presented study compares the AGIPD (pixel size of (200 jLim)^) to a possible apertured version 
of the detector and to a hypothetical system with (100 jitm)^ pixel size and investigates the influence of 
intensity fluctuations and incoherent noise on the quality of the acquired data. 

The intuitive conclusion that aperturing is not beneflcial as data is 'thrown away' was proven to be 
correct for low intensities. For intensities larger than approximately 1 photon per (100 jitm)^ aperturing was 
found to be beneflcial, as charge sharing effects were excluded by it. 

It was shown that for the investigated case (100 jitm)^ pixels produced signiflcantly better results than 
(200 |Ltm)^ pixels when the average intensity exceeded approximately 0.05 photons per (100 |Ltm)^. 

Although the systems were quite different in design they varied in the signal to noise ratio only by a factor 
of 2-3, and even less in the relative error of the extracted correlation constants. However the dependence on 
intensity showed distinctively different features for the different systems. 

Key words: AGIPD, simulation, XFEL, X-ray Photon Correlation Spectroscopy, aperturing, intensity 
fluctuations 



1. Introduction 

The European X-Ray Free Electron Laser (XFEL) [1 will provide ultra short, highly coherent X-ray 
pulses which will revolutionize scientiflc experiments in a variety of disciplines spanning physics, chemistry, 
materials science and biology. The European XFEL will provide pulse trains of up to 2700 pulses every 
222 ns (600 |Lts in total) followed by an idle time of 99.4 ms, resulting in a supercycle of 10 Hz and 27000 
pulses per second. 

One of the differences between the European XFEL and other free electron laser sources is the fast 
pulse repetition frequency of 4.5 MHz. Dedicated fast 2d detectors are being developed, one of which is the 
Adaptive Gain Integrating Pixel Detector (AGIPD) [2l|3ll5. The development is a collaboration between 
DESY, the University of Hamburg, the University of Bonn (all in Germany) and the Paul Scherrer Institute 
(PSI) in Switzerland. 
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12 The technique of X-ray Photon Correlation Spectroscopy (XPCS) provides an experimental method to 

13 probe dynamic properties of condensed matter [5 on nanometer length-scales. In XPCS experiments samples 

14 are illuminated by coherent X-ray light, which results in a grainy diffraction pattern (speckle pattern). As 

15 the speckle pattern is determined by the exact spatial arrangement of all the particles in the beam, any 

16 motion of the particles results in a fluctuation of the speckle pattern. The basic quantity determined in 

17 XPCS experiments is the intensity autocorrelation function of the speckles which measures tiny density 

18 fluctuations in the sample as a function of time and scattering vector Q. The high pulse repetition rate of 

19 the European XFEL, combined with the fast frame rate of the AGIPD, will allow investigations using the 

20 XPCS technique on the sub-microsecond timescale. 

21 Depending on the dynamic phenomena of interest and the associated time scales different technical 

22 realizations of photon correlation spectroscopy exist. For the very fast time scales in the regime of 10~^^ to 

23 10~^ seconds, the split and delay technique can be applied [610I8]. To study effects on the time scale of 

24 10~^ to 10^ seconds, which are relevant in soft matter or biological systems, sequential XPCS is performed 

- 0. 

26 Photon correlation experiments are limited by their signal to noise ratio (S/N or SNR) [9 . An important 

27 aspect of the SNR in XPCS experiments is the ability of the X-ray detector to resolve the speckle pattern. 

28 It can be shown that an optimal signal to noise ratio is achieved in situations where the speckle size equals 

29 the pixel size As XPCS seeks to measure dynamic sample properties, the X-ray pulses must illuminate 

30 the same sample spot several times. This implies that sample distortions (especially sample heating) due 

31 to interactions with the illuminating beam should be reduced as much as possible by using a large beam 

32 size. The speckle size is inversely proportional to the beam size and therefore speckle sizes will be small for 

33 XPCS experiments at XFEL sources. The pixel sizes of the fast 2D-detectors under development will be 

34 larger than the expected speckle sizes. This study investigates to what extent this mismatch will affect the 

35 SNR and the quality of the intensity correlation functions. Simulations of a prototypical XPCS experiment 

36 performed for different detector systems in otherwise identical situations are presented. The study is limited 

37 to small angle scattering experiments and shows the influence of the detector system on the SNR, which is 

38 commonly neglected. 

39 Furthermore, this study investigates a method of reducing the effective pixel size by aperturing the 

40 sensitive area of each pixel. Aperturing pixels produces a detector consisting of smaller pixels, that are 

41 isolated and separated from each other. Compared to adjacent small pixels, this eliminates charge sharing 

42 and makes it possible to sample a larger area with the same number of pixels. On the other hand the fraction 

43 of the detector area sensitive to photons is greatly reduced. 

44 This study is the first to investigate the performance of the AGIPD and the influence of certain design 

45 choices on the data quality of the experiment. A system level simulation study using an XPCS case for the 

46 DSSC[^ project can be found in literature [IJJ. 

47 2. Simulation of a prototypical XPCS experiment at the European XFEL 

48 All the simulations have been performed within the IDL framewor Inversion 7.1. The simulations consist 

49 of three independent stages modeling the real space system, the diffraction image generation and the detector 

50 response. The stages are explained in detail below. The simulations are followed by a data evaluation step, 

51 which is explained in section [3] 

52 The simulations use a photon energy of 12.4 keV (A = 1 A), which is the design energy of the European 

53 XFEL. Most XPCS experiments will be conducted using the Materials Imaging and Dynamics (MID) ex- 

54 perimental station, where the distance between sample and detector will be around 10 m. An illustration is 

55 shown in Figure [l] 
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Figure 1: Illustration of an exemplary XPCS experiment. Typical distances between sample and detector for experiments at 
the European XFEL are about 10 m. The grains on the diffraction patterns are speckles. The images used in this study differ 
from the shown image, as explained later in this manuscript. Image reproduced from the European XFEL TDR 

56 2.1. Real space system 

57 The simulated diffraction patterns were based on the evolution of a real space system. For ease of 

58 simulation the simulated real space system consisted of particles hopping on a two dimensional grid as 

59 illustrated in Figure |2] The particles did not interact with each other and could occupy the same grid 

60 position. Per simulation step (At) each particle could move by ±1 grid position with a probability of 1 % 

61 in each of the dimensions independently. Particles leaving the grid on one side were replaced by particles 

62 entering on the opposite side. In total 5000 particles were simulated on a square grid of 2000x2000 points. 

63 The initial distribution of particles was random. 

64 2.2. Generation of diffraction images 

65 For each simulation step the real space particle distribution was converted to an image. This real space 

66 image had the same size as the simulation grid and each image pixel was assigned the number of particles 

67 at the corresponding grid point. From this point onwards, all the calculations were done image-based. 

68 2.2.1. Speckle size, illumination function and Fourier transform 

69 The speckle size S is usually given by the following expression: 

S^'-§ (I) 



■"^The DEPFET Sensor with Signal Compression (DSSC) project and Large Pixel Detector (LPD) project are the other two 
projects developing a dedicated detector for experiments at the European XFEL. 

^IDL stands for Interactive Data Language and is distributed by ITT Visual Information Solutions. For more information 
visit [http : //www . ittvis . com/] 
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Figure 2: Illustration of the real space system. The particles did not interact with each other and could occupy the same grid 
position. Per simulation step each particle could move by ±1 grid position in each of the dimensions independently. Particles 
leaving the grid on one side were replaced by particles entering on the opposite side. In total 5000 particles were simulated on 
a square grid of 2000x2000 points. The initial distribution of particles was random. 



70 with A being the wavelength of the ihuminating hght, L the sample to detector distance and D the size of the 

71 illuminating light beam. Therefore an illumination function, which corresponds to the spatial distribution 

72 of the primary beam, is needed to define the speckle size in the detector plane. 

73 The complex wave field in the detector plane Ediff was defined as: 

Ed^ff=FT{Frs^Ful} (2) 

74 where F^s is the real space image and Fm the illumination function with FT denoting the discrete Fourier 

75 transform (FFT algorithm) and * indicating pixel by pixel multiplication. The illumination function was 

76 assumed to be a Gaussian function centered on the image center with a sigma corresponding to one eighth 

77 of the image size. 

78 To remove the structural artifacts of the square image the first 200 pixels in x and y direction were 

79 discarded. The removed area is shown in light gray in Figure [3j 

80 A Fourier transform of real values is symmetric around the Nyquist frequency, as indicated by Friedel's 

81 Law. As Ediff {Q) = E^-jj{—Q), areas symmetric to Figure [s] were removed from the complex wave field, 

82 leaving data from a region of 800x800 pixels for further processing, shown in medium and dark gray in 

83 Figure [3j The intensities from an exemplary wave field are shown in the left image of Figure [4] and exhibit 

84 no structure and form factor effects. 

85 The upper right image of Figure [4] shows a histogram of the intensities of the complex wave field with a 

86 bin size of one hundredth of the average intensity (/). The expected negative exponential distribution ^] 

87 for fully coherent illumination is clearly reproduced and the most probable intensity value is zero. The large 

88 variation in the histogram bins at large values of //(/) is due to the finite number of image pixels. 

89 The lower right image of Figure [4] shows the average normalized spatial autocorrelation function ACF 

90 of the intensities along the x direction averaged over the y direction defined by: 

ACF{AR} = ^ 2^ ^N^-i^j 7r^2 

91 where AR is the inter pixel distance, Nx and Ny are the number of pixels along the corresponding direction, 

92 //c,/ the intensity in the indicated pixel and {I)i the average intensity in the corresponding line. 

93 By calculating the average normalized spatial autocorrelation function the speckle shape and size was 

94 determined independently from the knowledge of the illumination function. The speckle shape corresponds 
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Figure 3: Regions of the diffraction pattern and their assignment in this study. Light gray: discarded pixels due to possible 
structure and form factors. Medium and dark gray: pixels used for all studies of (200 |Lim)^ detectors and small ROI investi- 
gations of (100 [xm)^ detectors. Dark gray: pixels used for studies of (100 |j.m)^ detectors and large ROI investigations. Pixels 
assigned to the exemplary Q bins of 300, 400, 500 and 600 units have been marked black. 

95 excellently to the Gaussian fit with a sigma of 1.25 pixel^ Exchanging x and y direction in the evaluation 

96 yielded identical results, indicating symmetric speckles. 

97 Speckle positions can be misaligned with the detector pixel grid. In order to simulate this, 2 by 2 pixels 

98 of the Fourier transform were assigned to each (100 jitm)^ pixel area. Correspondingly (200 jitm)^ area were 

99 assigned 4 by 4 pixels of the complex wave field. The sub-pixel position information gained by this procedure 

100 was used in the detector simulation step. 

101 2.2.2. Intensity normalization and fluctuation, image quantization and incoherent noise 

102 The diffraction pattern impinging on the detector consists of discrete photons. This discrete photon 

103 distribution was generated in the following way: 

^7 = ^1 {^)dRfluct I + ^noise (4) 

104 where S{X) is the image quantization procedure explained below, {I)d the desired average intensity, Rfiuct 

105 the intensity fluctuation factor explained below, I{Edi ff) = \Edi / / 1 the intensity distribution in the detector 

106 plane, {I{Ediff)) the average intensity in the detector plane and Fnoise the additional incoherent noise 

107 explained below. 

108 The fully coherent wave field is ergodicj^ In order to quantize the photon field (denoted by S{X))^ the 

109 photons were sampled according to the negative binomial distribution expected for its spatial distribution 
no [12], using the floating point pixel intensity as the expected value (/): 



^As the Fourier transform is oversampled and one pixel length corresponds to 50 fxm, this results in a sigma of 62.5 |Ltm and 
a full width at half maximum of 147.2 |j.m for the speckles in this simulation study. Assuming 10 m distance between sample 
and detector, this yields a FWHM of 6.8 jiim of the illuminating beam. 

"^In very simplified terms: An ergodic system shows the same average behavior when probed an infinite number of times 
sequentially as an infinite number of identical systems would show when probed once. 
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Figure 4: The left image shows the intensity distribution of one exemplary Fourier transform. No structure and form factor 
can be observed. The upper right image shows the histogram of the intensity distribution, binned with a bin size of 0.01 {/}. 
The lower right image shows the autocorrelation function along the x axis averaged over the y axis, thus revealing the average 
speckle shape. The speckle shape is excellently matched by a Gaussian fit with a sigma of 1.25 pixels. 
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^^^^-r(M)r(j + i)V V m) 



where / is a discrete number of photons, P{I) the probabihty of reahzing / photons and M the number of 
112 modes in the image. T{x) = (x — 1)! denotes the gamma function. Thus, for the fuhy coherent case (M = 1) 
investigated in this study, equation |5] simpUfies to: 



Pin = { ) (6) 



i + (j) Vi + (/) 

114 The intensity fluctuation factor Rfiuct was either set to unity (no fluctuations) or sampled randomly 

115 from a Gaussian distribution with a mean of 1 and an rms of 0.1. The factor does not simulate intensity 

116 fluctuations directly, but rather remaining uncertainties after accounting for intensity fluctuations. Intensity 

117 fluctuations can be accounted for in the determination of the autocorrelation function by weighting images, 

118 however the most probable value of zero photons cannot be scaled, hence creating uncertainties. Accounting 

119 for intensity fluctuations is particularly important at FEL sources, as the pulse to pulse variations can be 
quite large [13 . 

The incoherent noise Fnoise added after the quantization process was either zero or considered Poissonian 
with a probability of 1/400 per pixel of the Fourier transform. This translates to a noise photon probability 
of 1/100 per (100 jitm)^ detector area. Thus for an average intensity of 0.01 the SNR from the experiment 
itself is 1. As the noise of the detection process was simulated separately (see below) this completely random 
125 noise can be considered a high estimate of fluctuations of an otherwise constant background distribution. 
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126 2.3. Detector response 

127 In order to evaluate the effects of the investigated detector systems (detailed below), the same photon 

128 distributions were taken as input for the different detector systems. 

129 2.3.1. Simulated detector systems 

130 To evaluate the influence of different detector choices the following detector systems were simulated: 

131 • Ideal detector systems with (100 |Ltm)^ and (200 |Ltm)^ pixel size. 

132 • The Adaptive Gain Integrating Pixel Detector (AGIPD) featuring (200 jitm)^ pixel size. 

133 • A Modified AGIPD using Aperturing Techniques (MAAT), which was identical to the AGIPD, but 

134 completely insensitive to photons which are closer than 50 iitm to a pixel boundar}|^ which resulted 

135 in an effective pixel size of (100 jitm)^. 

136 • A Reduced Amplitude SEnsing System (RAMSES), which was a hypothetical detector which sacrificed 

137 most of the AGIPD functionality (e.g. the gain switching) but featured (100 |j,m)^ pixel size. 

138 2.3.2. The Adaptive Gain Integrating Pixel Detector (AGIPD) 

139 AGIPD is based on the hybrid pixel technology. A newly developed Application Specific Integrated 

140 Circuit (ASIC) will feature in each pixel a dynamic gain switching amplifier (to cope with the high dynamic 

141 range) and an analogue memory capable of storing approximately 300 images at the desired 4.5 MHz speed. 

142 The AGIPD will feature a pixel size of (200 jitm)^ and a silicon sensor with a thickness of 500 iitm. The 

143 area of (200 jitm)^ is needed to realize in-pixel electronics which fulfill the requirements for detectors at the 

144 European XFEL [2", "3^. The image data is read out and digitized in the gap between two pulse trains. 

145 As all the photons arrive at the detector in a period of time that is very short compared to the time 

146 needed for charge transport in the sensor material, the AGIPD finds the number of photons absorbed in 

147 each pixel by integrating the total signal. 

148 The AGIPD detector will be available for the first experiments at the European XFEL. The other 

149 systems studied here have been suggested as possible modifications for XPCS experiments. Preliminary 

150 studies showed that a system with reduced pixel size (RAMSES) would not be able to store as many frames 

151 as the AGIPD; maybe fewer than 100. 

152 2.3.3. Physical constraints precluding the realization of an ideal detector 

153 Due to the special pulse structure of the European XFEL, it is necessary to store the acquired images 

154 inside the pixel logic during the pulse train. A compromise has to be found between a large pixel area, so 

155 that many images can be stored, and a small pixel area for high spatial resolution. 

156 Additionally, effects of the limited stopping power of silicon and the charge transport in the sensitive 

157 volume come into play. The chosen sensor thickness of 500 |Lim provides a quantum efficiency exceeding 85% 

158 for 12 keV photons. After a photon is absorbed, a cloud of charge carriers is created, which expands before 

159 reaching the readout electrode and may cause charge sharing between pixels. The number of electron hole 

160 pairs created per absorption is normally distributed with a mean of (iVe,/i) = ^7/3.6 eV and a variance of 

161 ^ = F{Ne^h)^ with F ^ 0.1 being the Fano factor of silicon. 

162 The actual shape and size of this cloud is influenced by many parameter^ In this study, the collected 

163 charge is approximated as following a Gaussian distribution with a width depending of the square root of 

164 the distance to the electrode [TT. To minimize charge sharing between the pixels, the charge cloud size 



^For real applications this could e.g. be realized by post-processing a tungsten grid of 10 |Lim thickness on the entry window 
of the detector. 

non-exhaustive list includes: bias voltage, sensor doping, number of charge carriers, electrode geometry and exact initial 
position. 
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165 should be small compared to the pixel siz^ Given the material properties of a typical silicon detector of 

166 500 |LLm thickness and using recent determinations of the charge carrier mobilities p^, typical cloud sizes 

167 (full width at half maximum) are in the range of 20 jLLnj^ 

168 2.34. HORUS 

169 The detector response was calculated using the HORUS software described in ^ilTj. HORUS has already 

170 been successfully used to simulate the performance of the MedipixS chip [18l [19] and recently for AGIPD 

171 [2^. The entry window of the sensor was included as 3 iitm of insensitive silicon equivalent material for all 

172 systems, which is about twice the thickness foreseen for AGIPD at the moment. For all systems, the total 

173 noise was assumed to have an equivalent noise charge of approximately 300 electrons, which is about double 

174 the value determined from ASIC simulations [3 and roughly the value measured on a test assembly using a 

175 microstrip sensor. 

176 To increase the simulation speed several approximations were used: 

177 • Any parallax effects were ignored. This is a reasonable approximation, as the sample to detector 

178 distance for XPCS experiments at the European XFEL will probably be 10 m or more. 

179 • Any effects originating from sensor edges and module gaps were ignored. Although in the final detector 

180 system these effects cannot be avoided, for XPCS data evaluation it is possible to mask areas exhibiting 

181 these effects and exclude them from the data analysis. 

182 • Effects produced by high instantaneous charge densities (so called plasma effects) are not simulated 

183 by HORUS, and thus were not included in this study. Neglecting these effects is a reasonable approx- 

184 imation, as the intensities in XPCS experiments are usually low and plasma effects are expected to 

185 show only for intensities larger than 250 photons per pulse per (100 jitm)^ area j2T| [22| [23]. 

186 The ideal detectors were realized by taking the total number of photons impinging on each pixel. Thus 

187 there is no detector noise and the quantum efficiency of the ideal systems is 1.0 (in contrast to ~ 0.89 for 

188 the non-ideal systems). 

189 The response of the MA AT was calculated in an identical way as for the AGIPD, but photons hitting 

190 the insensitive regions were discarded in F^. 

191 The RAMSES was modeled using the AGIPD parameters, except that the pixel size in the simulations 

192 was (100 |Ltm)^, quadrupling the number of pixels in the simulation (as the same Q space is covered). 

193 2.3.5. Region of interest 

194 There will be experimental cases where the entire region of interest can be sampled with the detector 

195 area. For these cases the detector performance is limited by the pixel density, and similar total areas covered 

196 by (100 |Ltm)^ and (200 |Ltm)^ pixels had to be compared. These cases are labeled as cases with small Region 

197 Of Interest (sROI) in the rest of this study. An illustration of this is shown in Figure [s] 

198 For cases where the potential region of interest exceeds the detector area, (100 M-m)^ and (200 M-m)^ 

199 systems with an identical number of pixels had to be compared. The Q region for the simulated small pixel 

200 detectors was limited to the same number of pixels that was used for the large pixel detectors, shown as the 

201 dark gray area in Figure |3l These cases are labeled as cases with large Region Of Interest (IROI). 

202 It should be noted that, translated to real physical systems, the Q vector units in the small and large 

203 ROI cases would therefore have different meanings. 



^For XPCS applications charge sharing should be avoided if possible, as the signal of adjacent pixels would be correlated, 
increasing the effective pixel size. It should be noted that for very low intensities (low pixel occupancies) charge sharing is 
beneficial, as event-by-event correlations can be performed [15| . 

^200 V sensor bias results in drift times td around 30 ns, assuming the whole thickness as drift distance. The lateral cloud 
size can then be approximated as FWHMdoud ~ 2.35v^2Dt^ and D = fikT/e with fi,k,T,e being hole mobility, Boltzmann 
constant, temperature and electron charge, respectively. The FWHM value of 20 fxm is the result of assuming (100) crystal 
orientation and a temperature of 300 K. 
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Figure 5: Illustration emphasizing the difference between small and large areas of interest. The actual pixels simulated and 
their assignment to the regions of interest is explained in detail in Figure [3] 



204 3. Data evaluation 

205 The simulation was performed multiple times for each intensity, simulating the data acquisition of five 

206 independently pulse trains {Ntr = 5), thus the overall data set consisted of 5 x 300 frames {Nf = 300) per 

207 detector type per intensity point. 

208 3.1. The intensity autocorrelation function 

209 The intensity autocorrelation function (^2 function) is both a measure of the correlation time and the 
contrast. For ideal circumstances, which are not realized here, the ^2 function would be a function mono- 
tonically decreasing from a value of 2 for kAt ^ to a value of 1 for kAt 00. The functional form 
is determined by the sample dynamics and the correlation time constant can be extracted from this. The 
scaling is a measure of the optical contrast. 

For the data evaluation the ^2 function was calculated individually for each pixel p according to the 
following formulap^ 



, 1 1 ^'^^ Ip^tr{iAt)Ip^tr{{i^k)At) 

g^ip.kAt)-^ ^ 2^ 



tr 



tr=0 



i=0 



where the average intensity of pixel p {{Ip)) is calculated as: 



Nf-1 



{Ip) 



tr=0 f i=0 



(7) 



(8) 



217 and At is the time difference between frames and Ip^tr {iAt) denotes the intensity in pixel p in frame i of 

218 pulse train tr. This procedure yielded a ^2 function for each pixel of each investigated average intensity and 

219 detector type. 



^In the simulations each train used a new random number sequence. 

^"-"The shown equations use the commonly employed normalization scheme. The data was also analyzed using the symmetric 
normalization scheme discussed in [24], however the results were worse then for the commonly employed normalization scheme 
and are therefore not shown here. 
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Figure 6: Number of pixels as a function of Q value for the simulated detector systems. The Q range between the dashed lines 
has been used for Q averaging, as explained in the text. The dotted line is the Q vector shown in Figures |7| [s] |13| and |14| 
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Figure 7: Resulting g2 functions at a fixed Q of 500 units as a function of average intensity. Results were obtained for the 
ideal (100 ixm)-^ detector for small ROIs. The results include photon statistics noise and noise from background photons and 
intensity fluctuations. The apparent increase of the g2 at very large lag times is an artifact induced by the flnite number of 
images. The noise is especially prominent at the lowest intensities. 



220 Especially at very low intensities, there were pixels which did not detect any photon in the whole series. 

221 To avoid division by zero (as (Ip) = 0) the g2 function for these pixels was set to zero. 

222 3.2. Q space binning and azimuthal average 

223 Azimuthal averages (ensemble averages) were calculated by binning pixels that fulfilled the condition 

224 Q — AQ < Q{p) < Q + AQ, where Q{p) was the absolute value of the Q vector for pixel p and 2AQ = 10 

225 was the bin size. Exemplary Q space bins are shown in Figure |3j the number of binned pixels as a function 

226 of Q value is shown in Figure [gJ The ensemble averaged g2{Q^ kAt) function is the average of all g2{p, kAt) 

227 functions in the corresponding Q space bin. 

228 As an example, the resulting ^2(Q, kAt) functions at a fixed Q value of 500 units as a function of average 

229 intensity are shown for the ideal (100 jitm)^ detector in Figure and for AGIPD and MA AT in Figure p] 
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Figure 8: Resulting §2 functions at a fixed Q of 500 units as a function of average intensity. The upper image shows the results 
for the AGIPD, the lower image shows the results for the MAAT. The results include photon statistics noise and noise from 
background photons and intensity fluctuations. The apparent increase of the g2 at very large lag times is an artifact induced 
by the flnite number of images. The noise is especially prominent at the lowest intensities. 



230 Figure [7| shows the expected exponential decay of the autocorrelation function with increasing lag time. 

231 The maximum value for low lag times is approximately 1.7, thus below the theoretical maximum value of 2, 

232 as the experimental contrast is smaller than one. For lag times beyond 100 the function falls below 1.0 and 

233 increases again towards larger lag times. This behavior is commonly observed when evaluating the intensity 

234 autocorrelation function on a limited set of frames and would disappear if more frames were imaged in 

235 sequence. The apparent contrast decreases as the intensity decreases. At the lowest intensities the images 

236 become uncorrelated as they are dominated by noise, correspondingly the g2 function appears flat. 

237 Figure [S] shows the exponential decay of the autocorrelation function for the AGIPD and MAAT. The 

238 features seen for the ideal system are observed here as well. The maximum value for AGIPD and MAAT at 

239 low lag times is approximately 1.3 and 1.7, thus MAAT produces a similar contrast to the ideal (100 |Ltm)^ 

240 detector. When compared to the g2 functions shown in Figure [7| the AGIPD curves appear smoother, as 

241 the larger pixel size averages over a larger area, while the MAAT curves appear less smooth, as the the g2 

242 function is calculated from less detected photons. This is especially prominent at the lowest intensities. 

243 3.3. Fit to the g2 function 

244 Each ensemble averaged g2{Q-, kAt) function was fitted with the following function to evaluate the per- 

245 formance of the different detector systems: 
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Figure 9: Average scaling parameter for all systems. The left shows the results for photon statistics noise only, the right image 
shows the results including photon statistics noise and noise from background photons and intensity fluctuations. The dashed 
vertical line indicates the average intensity of Fnoise- 



52 {Q, kAt) = S{Q) {C{Q) exp(-2r(Q)T) + 1} (9) 

246 where r = kAt is the lag time, S{Q) a scahng parameter, C{Q) the optical contrast and T{Q) the inverse 

247 of the correlation time Tc. 

248 For the dispersion relation T{Q) it was assumed that T{Q) = CdispQ^ and likewise Tc = {CdispQ)~^ • 

249 The scaling parameter S{Q) was necessary to compensate the effect of pixels which detect zero photons 

250 during the bunch train. It can be calculated analytically as 1 — {^N^^\^^q/Nq) ^ where N^j^^^q/Nq is the 

251 fraction of pixels in a given Q bin which have zero intensity in all frames of a given pulse train. 



252 4. Results 

253 The influence of the intensity fluctuation factor Rfiuct in Equation |4] on the results was found to be 

254 negligible. Thus two sets of results are compared here: the results with photon statistics noise alone 

255 {Rfiuct = 1.0 and Fnoise = 0), and the results with photon statistics noise and noise from background 

256 photons and intensity fluctuations. 

257 Using the Fourier transform to generate diffraction patterns implies an infinite coherence length, both 

258 longitudinal and transverse, thus the commonly observed drop in contrast at high Q values is neglected. As 

259 a consequence there is little change in contrast and scaling factor as a function of Q value. Thus the fit 

260 results were averaged for Q values between 310 and 520 units. 

261 The Q range between 310 and 520 units has been chosen as it excludes Q vectors with low statistics 

262 towards the edge of the sensitive area for large regions of interest. The (200 jitm)^ pixel detectors and small 

263 ROI cases have been restricted to the same Q range to have a comparable data set, although the errors for 

264 these systems could have been reduced using the additional information available at other Q values. 



4.1. Scaling parameter 

Generally the scaling parameter shown in Figure [9] behaves as expected and explained in section |3.3 



267 The average scaling parameter for the (200 jam) systems is always larger than for the (100 jam) system as 

268 the probability to measure photons increases with the sensitive area. 

269 When summing together all the images of a pulse train, the scaling factor can be interpreted as the 

270 fraction of non-zero pixels in the summed image. This can be clearly observed as without noise, the scaling 

271 factor is decreasing at low intensities. 
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Figure 10: Average optical contrast for all detector systems as a function of intensity. The left image shows the results for 
photon statistics noise only, the right image shows the results including photon statistics noise and noise from background 
photons and intensity fluctuations. The dashed vertical line indicates the average intensity of Fnoise- 



272 For the summed image the expected number of noise photons N^^noise can be analyticahy calculated to 

273 be N^^noise = / Pnoise = 3.0 photous per (100 jLim)^. Since these photons follow Possionian statistics, the 

274 probability to have zero photons in a pixel of the summed image when the expected value is 3 is ^3(0) ~ 0.05. 

275 When noise is present an average scaling factor of 1 — ^3(0) ~ 0.95 is expected for (100 jitm)^ systems, and 

276 this is observed for low intensities. 



4.2. Optical contrast 



Figure [To] shows the Q averaged results for the optical contrast. When Fnoise is included in the simulations 
the optical contrast decreases once {I) /Inoise ^ 100 ((/) ^ 1), while it stays almost constant without noise. 
For Fnoise = and very low intensities the contrast apparently increases due to the biasing effect of the 
scaling factor. However this additional contrast did not improve the data quality, as the scaling factor 
decreased strongly. 



Figure 10 shows that, except for MAAT, the real systems show a lower contrast than the ideal ones. 



independent of the noise. The effect is small for the AGIPD, as the vast majority of photons do not 



undergo charge sharing, but the observed small difference between the green and yellow lines in Figure 10 
is statistically significant. At low intensities the difference vanishes. This reduced contrast was caused by 
charge sharing events, which correlate the pixel values. Charge sharing is excluded for MAAT by its design. 

The fact that the results for MAAT reproduce the results for an ideal (100 jitm)^ system, indicates that 
the detector noise and quantum efficiency do not influence the data quality. 

For the ideal systems the contrast was given by C((/)) ^ Cgeom{^ — Inoise/ (^)), where C is the determined 
contrast and Cgeom the contrast expected from geometric estimations: 

Cgeom ~ ^2" (10) 

with P being the linear pixel size and S the full width at half maximum of the speckles. 

The geometric estimate and the results of the ideal systems at the highest intensity are in good agreement 
(0.68 and 0.72 or 0.35 and 0.36 for (100 |Lim)^ or (200 |Ltm)^ systems, respectively). 



4.3. Correlation time 

shows the inverse correlation time T{Q) for the ideal (100 jitm)^ pixel detector and small ROIs 



Figure 
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as a function of Q vector and intensity. 
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Figure 11: Inverse correlation time r(Q) for the ideal (100 |j.m)^ detector and small ROIs as a function of Q vector and intensity 
for photon statistics noise only. The results for the other detector systems look similar, but show larger variations. 
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Figure 12: Dispersion constant C^isp for all detector systems as a function of average intensity. The left image shows the results 
for photon statistics noise only, the right image shows the results including photon statistics noise and noise from background 
photons and intensity fluctuations. The dashed vertical line indicates the average intensity of Fnoise- 
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Overall, the correlation time was observed to be approximately inversely proportional to the Q vector 
(r(Q) = CdispQ-, likewise Tc = (CdispQ)~^)- The Q averaged results for the dispersion constant Cdisp are 



displayed in Figure [12] for all simulated detector systems. 

For purely diffusive systems T{Q) should be proportional to Q^, but the investigated jump diffusion 
system did not show such a behavior. Literature shows [25l [26] that a diffusive behavior is only expected 
for low Q vectors, where the length scale is large compared to the jump size. 

As the low Q vectors are not included in the data evaluation (detailed in 2.2.2), a dependence is 
neither observed nor expected. 

Figure 12 shows a slightly different dispersion constant Cdisp foi" each system. This is due to the different 



Q space coverage of the systems and different averaging properties of the different pixel sizes and the crude 
approximation of an inversely proportional dispersion relation, which is purely empiric. However the stability 
of the results as a function of intensity show that this is merely a systematic error, which is identical for all 
systems of the same type, so it will not influence possible conclusions. 



When evaluating the influence of the average intensity in Figure 12 one notices instabilities for the 

14 



all 3 noise sources 




Figure 13: Signal to noise ratio for all detector systems as a function of average intensity including photon statistics noise and 
noise from background photons and intensity fluctuations. The curves shown have been determined for a Q value of 500, other 
Q values produce similar results. The dashed vertical line indicates the average intensity of Fnoise- 



results at average intensities < 0.01. At these intensities the fit results are of low statistical significance. 
The extracted dispersion constant Cdisp was observed to be more sensitive to the low significance than the 
contrast. A reason for this might be the additional data evaluation step necessary to extract Cdisp- 

When evaluating the influence of the added incoherent noise by comparing the left to the right image 



in Figure 12 at average intensities > 0.02, one notices very little effect. The reason of this is that the 
correlation time is independent of the contrast, as long as the contrast is high enough to reliably determine 
the correlation time. 

5. Signal to noise ratio 

The signal to noise ratio is commonly used as a figure of merit for XPCS experiments. Adapted to the 
situation of this study the analytic expression for the SNR, derived from the mean value of the g2 function 
and its error (err(^2)J3 is: 

sNRm = oc {i,)cm^ (11) 

with (Ip) being the average intensity in a pixel area, C{{I)) the optical contrast at this intensity and N the 
number of detector pixels. 



Figure 13 shows the signal to noise ratio determined by dividing the first value {kAt = 1) of the g2 



function minus 1 by its error determined from the variance of all g2 functions at this lag time in the same 



-'^-'^The error of the mean X using Ns samples of a statistical population X is usually deflned as err(X) = ^Jvars{X) / Ns^ 
where vars{X) = {^i ~ ^)'^ ^^e variance of the Ng samples of X. 
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Figure 14: Signal to noise ratio for selected detector systems in the low intensity regime. The left image shows the results 
for photon statistics noise only, the right image shows the results including photon statistics noise and noise from background 
photons and intensity fluctuations. The dashed vertical line indicates the average intensity of Fnoise- The analytical SNR for 
(200 |xm)^ systems is almost identical to the analytical SNR for (100 |xm)^ systems and small regions of interest; their lines 
coincide almost completely. 



Q space bin. This is the definition of the signal to noise ratio in the hmit of smah lag times as indicated in 
Equation [iT] Figure 13 shows the signal to noise ratio at a fixed Q value of 500 (indicated by the dotted 
line in Figure [6|. Evaluation of the signal to noise ratio at different Q values showed identical results. 

The main features of Figure 13 are the increase in SNR with intensity for low intensities and the saturation 
of the SNR for high intensities. The saturation value is reached around average intensities of 1 photon per 
(100 |Ltm)^. Additionally Figure 13 shows multiple separate sets of lines, one set for (100 |Ltm)^ with small 
ROI, one set for large ROI and one set for (200 jitm)^ systems. For low intensities the signal to noise ratio 
of the (200 jLim)^ systems is approximately equal to the SNR of (100 jitm)^ systems with small ROI and 
approximately twice the SNR of the (100 jitm)^ systems with large ROI. Towards saturation all (100 jitm)^ 
systems show larger SNRs than the (200 jitm)^ systems, although only a little for large ROI and MAAT. 
The ratio of the SNR between large and small ROIs corresponds the square root of the ratio of the number 
of pixels. Ideal systems produce an SNR which is a little higher than that of the real systems due to their 
perfect quantum efficiency. 

Equation [TT] was derived under the assumption of low count rates and small contrast in the limit of small 
lag times ([9l[T0] and references therein). A more elaborate derivation, also for large lag times, can be found 
24 . According to Equation [ill the SNR increases with intensity until the assumptions for its derivation 



m [:^4 . According to Equation 
are no longer valid. 



344 Using Equation 11 to compare the SNRs of two systems with linear pixel sizes Pi and P2 yields: 



SNR^{{I)) (Ji) Ci((J)) 



(12) 



SNR^m ih) C2({/)) 
Assuming that the intensity in the pixel is proportional to its area and the intensity dependent terms of 



346 the contrast cancel, one can rewrite Equation n2\ using Equation 10 



as: 



SNRi 

SNR2 ~~ P| 52 + p2 




(13) 



Equation 



13 



yields values from y/N^/N^ for Pi,P2 > ^ to (P1/P2) for Pi,P2 <C S. For the 
systems investigated here Equation [13] yields the results observed for low intensities, namely 

349 SNRioo^rn,sROi/SNR2oof^m ^ 1-0, S N Riqq ^^n ,1 RO I / S N R200 ^ 0.5 and 

350 SNRioo^rn,sROl/SNRioo^rn,lROI ^ 2.0 
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Figure 15: Relative error E of the dispersion constant C^isp for all detector systems as a function of average intensity. The left 
image shows the results for photon statistics noise only, the right image shows the results including photon statistics noise and 
noise from background photons and intensity fluctuations. The dashed vertical line indicates the average intensity of Fnoise- 




0.10 



Figure 16: Relative error E of the dispersion constant Cdisp for selected detector systems in the low intensity regime. The left 
image shows the results for photon statistics noise only, the right image shows the results including photon statistics noise and 
noise from background photons and intensity fluctuations. The dashed vertical line indicates the average intensity of -f^noise- 



Figure [M] shows the results for selected detector systems in the low intensity regime compared to the 
results from Equation [TT] It can be seen that Equation [TT] is a well suited approximation of the SNR at low 
intensities. Furthermore it can be seen, that for intensities below approximately 0.05 photons per (100 jitm)^ 
the signal to noise ratio of the AGIPD is larger than estimated from Equation [TT] 



6. Relative error E of the dispersion constant Cdisp 

The focus here is on a slightly different aspect of the error determination in XPCS experiments. As the 
main physical quantity of interest are the time constants connected with the dynamics the relative error in 
the time constants given by the ratio E — var{Cdisp) / Cdisp was investigated. This quantity was derived 
by fitting the g2 functions with the exponential correlation functions (Eq. [9|. Thus E represents a different 



figure of merit than calculated by Equation 11 While for the SNR higher values indicate better quality, 
lower values of E indicate better values. 

Both images of Figure 15 show a similar behavior in dropping from very large values at low intensities 
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to a saturation value at high intensities. As already observed for the correlation time, the errors are very 
large for intensities < 0.01, indicating problems with the proper determination of the resulting quantities. 

Comparing the left to the right image reveals the influence of the additional noise. While the data 
without additional noise is still usable at very low intensities the error increases significantly when the 
average intensity of the signal becomes comparable to the average noise intensity. 

For average intensities above approximately 0.1 the results are independent of the noise. From this point 
on the error contributions of the photon statistics and all other noise sources are small compared to the 
systematic errors. The systematic errors are introduced by artifacts in the g2 functions causing a small, but 
systematic deviation of the ^2 functions from the function used for the fit. The artifacts in the ^2 functions 
are introduced as the correlation time is a significant fraction of the total sampling time, which is limited 
by the total number of framej^ 

The saturation values are qualitatively different from the saturation values of the signal to noise ratio. 
The relative error saturates at lower intensities than the signal to noise ratio and to qualitatively different 
values, which result from the different numbers of pixels for the different detector systems, but do not show 
a \/N behavior. The evaluated data is restricted to the Q range of 310 - 520 (shown in Figure [6|, so the 
results were biased in favor of the (100 jitm)^ detector systems, as more pixels were used for the evaluation 
than for the (200 |Ltm)^ detectors. 

Figure [16] shows the same data as Figure [15) but for low intensities on a linear scale. 



Without additional noise {Fnoise — O7 l^ft image of Figure 16^ the difference in E between the systems 
is small, especially compared to the situation including additional noise (right image of Figure 16). 

It is observed that for (/) ^0.05 the relative error E of AGIPD is lower than the relative error determined 
for RAMSES and large regions of interest independently of the noise. While the significance of this result 
is debatable for Fnoise — 0, the result is statistically significant when additional noise is present. It should 
be noted that the higher signal to noise ratio of AGIPD does not translate directly to a lower relative error. 

7. Summary and discussion 

A simulation tool using the detector simulation program HORUS was set up in order to simulate XPCS 
experiments. The numeric simulations took the whole chain into account: modeling a simple real space 
system, generating diffraction images, finding the detector response and evaluating the data. It was shown 
that the simulated diffraction patterns show the counting statistics of fully coherent speckle patterns and 
have a well defined speckle size. 

The well known intensity autocorrelation function was used as a method for data evaluation and the 
most relevant parameters were extracted by a simple exponential fit to the ^2 function. The simulations 
were performed for a set of different detector systems and a set of different noise sources. 

A decrease of optical contrast was observed for the non ideal systems and attributed to charge sharing 
effects. 



The simulated intensity fiuctuations (see section 2.2.2) had a very small effect on the data evaluation 



which proved the robustness of the intensity autocorrelation technique and its applicability for experiments 
at free electron laser sources. 

Contrary to common expectation (see Equation [IT]) the signal to noise ratio was found to saturate above 
intensities larger than approximately 1 photon per (100 jitm)^. Deviations from the analytical expression 
were already observed more than one order of magnitude below this intensity. 

The signal to noise ratio observed for the AGIPD showed an unexpected transition behavior. At low 
intensities ((/) < 0.1) it is similar the signal to noise ratio of RAMSES and small regions of interest (limited 
by pixel density). At large intensities ((/) > 1.0) it is similar to the signal to noise ratio of RAMSES and 
large regions of interest (limited by total sensitive area). 



■"^^A possible remedy for this might be the use of alternate (e.g. logarithmic) sampling schemes, which will be investigated in 
a future study. 
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MAAT 



RAMSES IROI 



RAMSES sROI 



key 
features 



effective (100 |Ltm)^ pixel 
size due to aperturing 



(100 jLim)^ system with the 
same number of pixels 



(100 |LLm)^ system with the 
same area 



contrast 



larger, reproducing ideal 
(100 jLim)^ system 



larger, but lower than ideal 
(100 |LLm)^ system for 
(/) > 0.1 due to charge 
sharing 



like RAMSES IROI 



signal to 
noise ratio 



lower for (/) < 1.0, slightly 
larger otherwise 



like MAAT 



lower for (/) < 0.05, larger 
otherwise 



relative 
error E 



larger for (/) ^0.1, about 
the same otherwise 



larger for (/) < 0.05, 
although hardly significant 
without additional noise, 
lower at larger (/) 



lower except for the two 
lowest intensities, although 
hardly significant there 



Table 1: Performance results presented in this study in comparison to the AGIPD system. Intensities are given as photons per 
(100 |j.m)^ area. 



408 The error in the relevant determined physical parameters was shown to deviate in behavior from the signal 

409 to noise ratio. The relative error saturated as intensity increased, but at intensities above approximately 

410 0.1 photon per (100 jitm)^. This is about one order of magnitude lower than the saturation of the signal to 

411 noise ratio. 

412 A comparison of the performance of the detector systems is summarized in Table [l] The results for the 

413 signal to noise ratio of MAAT were partially anticipated in [27 . Although the systems are quite different 

414 in design they vary in the signal to noise ratio only by a factor of 2-3, and even less in the relative error. 

415 However the dependence on intensity shows distinctively different features for the different systems. 

416 The amount of simulated data (5 pulse trains) limited the precision of the results at the lowest intensities. 

417 Additional constraints were imposed by the limited number of acquired frames per pulse train. The system- 

418 atic infiuences of this proved to be a restricting factor. The finding that the number of frames infiuences the 

419 precision is in good agreement with theoretical predictions (see. e.g. Equation (18) and (19) in [28]) and is 

420 based on the fact that the g2 function measures pairs of intensities (i.e. correlations). Obviously, the more 

421 number of pairs (i.e. frames) are available in an experiment the better the statistics of the measurement. 

422 The simulation was based on the assumption that all detector systems provide the same number of stored 

423 frames in a pulse train. In reality a detector system with smaller pixels may lead to a reduced number of 

424 stored images. This is affecting both the time scale that can be covered in an XPCS experiment and the 

425 SNR in the limit of low intensities. 

426 In the simulations the sample dynamics were completely independent of the illumination process. In 

427 experiments some of the incoming photons will always be absorbed by the sample and its suspending 

428 medium, giving rise to a temperature increase and possible changes in the sample. These effects were 

429 completely neglected in the simulations presented here, but might ultimately be the dominating effect for 

430 experiments at the European XFEL. 

431 Furthermore, in experimental situations many more complications may occur, e.g. partial coherence and 

432 anisotropic speckle sizes, which have been neglected here. The impact of those effects is beyond the scope 

433 of this study, but might infiuence the results for experiments at the European XFEL. 

434 XPCS experiments require the positional jitter of the illuminating beam to be small compared to the 

435 beam size. It has been recently demonstrated at the Linear Coherent Light Source (LCLS) at Stanford, 

436 USA, that a jitter between 4% and 20% of the beam size can be achieved [29 , thus giving confidence that 

437 postional jitter will not seriously infiuence XPCS experiments at the European XFEL. 
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438 8. Conclusions 



439 The intuitive conclusion that aperturing is not beneficial as data is 'thrown away' was proven to be 

440 correct for low intensities (see Table [T]). For intensities larger than approximately 1 photon per (100 jLtm)^ 

441 aperturing was found to be beneficial, as charge sharing effects were excluded by it. 

442 It was shown that for the investigated case (100 jLtm)^ pixels produced significantly lower relative errors 

443 on the extracted correlation times than (200 jitm)^ pixels when the average intensity exceeded approximately 

444 0.05 photons per (100 |Ltm)^. 

445 The signal to noise ratio of the investigated (100 jitm)^ pixel systems exceeded the signal to noise ratio of 

446 the AGIPD when a) the average intensity exceeded approximately 0.05 photons per (100 |Ltm)^ and the same 

447 area is covered by both systems or b) average intensity exceeded approximately 1 photons per (100 |Ltm)^ 

448 and both pixels have the same number of pixels. 

449 It was found that under certain conditions the maximum number of frames limited the precision. Thus 

450 for XPCS experiments at the European XFEL maximizing the number of stored frames per pulse train 

451 might be as important as minimizing the pixel size, which are conflicting requirements. The use of alternate 

452 (e.g. logarithmic) sampling schemes to compensate for this limitation will be investigated in a future study. 
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